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ffOCKDim  MOB  BUMUNOT  mK>D 


I . INTRODUCTION 


Schlichting^  presented  a solution  for  a laminar  (or  turbulent)  gas 
jet  exiting  from  a small  circular  orifice  and  mixing  with  the  surround- 
ing fluid.  The  pressure  along  and  across  the  jet  is  assumed  constant. 

The  flow  is  axisymmetric  and  steady.  Boundary  layer  theory  is  applicable 
to  this  problem  and  the  velocity  profiles  are  considered  to  be  similar. 

Under  the  similarity  transform,  the  governing  equation  is^ 


FF*  _ F^  _ FF"  _ d_  /pit  _ F^\ 
2 n n dn  ' n / 


where 


prime  (')  denotes  differentiation  with  respect 
to  r\, 

ip  is  the  stream  function, 

V is  the  kinematic  viscosity, 

X is  the  axial  coordinate,  and 
y is  the  radial  coordinate. 

Equation  (1)  is  a third  order,  ordinary,  nonlinear  differential  equation. 
The  boundary  conditions  are^  at  n = 0,  F = 0 and  F'  =0. 

An  exact  form  of  Schlichting's  solution  is  obtained  for  the  given 
boundary  conditions  and  an  exact  solution  is  obtained  for  a more 
general  class  of  problems  where  the  boundary  conditions  are  not  restricted 
to  those  specified  for  a freely  expanding  jet. 

II.  EXACT  SOLUTION  OF  THE  CIRCULAR  GAS  JET  EQUATION 


Equation  (1)  is  integrable  since 


_d_  fFF'\  _ _ FF' ' f'^  ^ FF^ 

dn  \ n / ~ n n _2  ’ 


and  direct  integration  of  Equation  (1)  yields 


Sahliahting,  Boundapu  Lauer  Theorit,  MaGraw-Hill  Book  Co.,  6th  ed., 
N.y.,  1968,  p.  219-221. 


(2) 


FF'  - F”  i-  + r 
n n 1 


ff'  = f'  - nF", 


for  the  given  boundary  conditions.  Equation  (3)  is  a second  order, 
nonlinear,  ordinary  differential  equation. 

Now,  Equation  (3)  may  be  integrated  since 


^ (nF')  - 2F’  = nF"  - f’. 


= FF  , 


or  Equation  (3)  becomes 


2nF'  = 2C2  + 4F  - 


It  is  noteworthy  that  a solution  to  Equation  (3)  is  given  by  Kamke^ 
by  transforming  the  independent  variable.  Since  direct  integration  is 
possible,  such  transformations  are  unnecessary. 


It  is  possible  to  separate  the  variables  in  Equation  (4)  or 
2dF  _ dn  


2C2  + 4F  - F 


2 n 


and  integration  yields. 


1 ^ Si.-!  2 - F 
_ tanli  -^3^  = tn  n + C^, 


where 


X = ± 


+ 4 


Then  the  function  F becomes 

2 - F 
2x 


tanh  I Jin  (nC^) 


CnCj)"’'  - (nCj)"' 
(nCj)"''  + (nCj)’'  ' 


^E.  Kamke,  Differentialgleiahungen  Loeungmethoden  und  Loaungen,  Chelsea 
Publishing  Co.,  1948,  p.  562. 
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or 


2x  (nC  - 2x  + 2 (nC  ♦ 2 

F = 

1 + 

is  an  exact  solution  of  Equation  (3). 

For  the  freely  expanding  jet,  atn=0;F=0  and  F'  = 0 or  from 
Equation  (4),  C2=0orx=±l.*  Then, 

4 (nC 


1 + (nC3)^ 

represents  an  exact  solution  to  Equation  (3),  where  is  an  arbitrary 
constant.  This  solution  satisfies  the  boundary  condition  that  F = 0 
and  f'  = 0 at  n = 0,  for  any  C3.  Note  that  for  C3  = where  y is  a 
constant,  and  for  C = Yb,  F transforms  to  Schlichting's^  solution  or 


However,  any  multiple  of  C3  will  yield  a similar  expression  for  F which 
satisfies  the  posed  boundary  condition.  The  general  solution  is  given 
by 

4 (nC,)^ 

F = 2 > 

1 + (bCj)^ 

and  the  second  constant  of  integration,  C3,  may  be  determined  directly, 
from  the  given  value  of  the  momentum  in  this  case. 

III.  EXACT  SOLUTION  OF  THE  GENERAL  EQUATION 

The  general  expression.  Equation  (2)  may  also  be  solved  exactly 
when  it  is  not  required  that  Cj  be  zero.  Equation  (2)  can  be  written 
as 

FF'  + nF"  - f'  - n Cj  = 0, 
and  integrated  to  yield 

F^  + 2nF'  - 4F  = + 2C2.  (5) 

* x = +1  and  X = -1  yield  solutions  of  exactly  the  same  form. 


I 
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r 


2 2 

Now  F - 4F  = (F-2)  - 4,  and  Equation  (5)  becomes 

(F-2)^  + 2riF'  = + 4.  (6) 

The  left  hand  side  (LHS)  of  Equation  (6)  is  simplified  by  letting 
“(n)  = F(ti)  - 2,  and  Equation  (6)  becomes 

+ 2na'  = + 20^  + 4.  (7) 

The  term  n is  removed  from  tl^e  LHS  of  (7)  by  transforming  the  independent 
variable  n such  that  a(n)  = a(C)  and  5 = iln  |n|,  or 


^ + 2a'(5)  = 2C2  + 4 + Cj  e^^. 


(8) 


The  constant  coefficients  are  rearranged  by  letting  y(C)  = ^ , 
then  Equation  (8)  becomes. 


e^^  C,  C 

y'  • y'  '■ 


(9) 


Equation  (9)  is  of  the  Riccati  type  and  may  be  transformed  to  a linear 
equation.  This  is  accomplished  by  the  Riccati  transformation 


and  Equation  (9)  becomes 


u'(C)  = y(0  u(C), 


'e^'c  C 

u [—  * T * ‘ ■ 


(10) 


Now,  the  inverse  transformation  on  the  independent  variable  C,  where 

/c7 


u(C)  = u(z) 


v: 


z = w/ 


yields. 


2 / 2 ^2  , , 

Z U +ZU  -Iz  +-^+l)u=0. 


(11) 


2 S 

Equation  (11)  is  the  modified  Bessel  equation  for  v = -^  + 1 . The 
solutions  are  I^(z)  and  K^(z)  which  are  modified  Bessel '^functions  of 


the  first  and  second  kind  respectively^,  or 


u = A 1^(2)  + B K^(z).  (12) 

where  A and  B are  arbitrary  constants  and  1^  and  I_^  are  linearly 
dependent  if  v is  an  integer.  Also  I^Cz)  and  K^(z)  are  real  and  posi- 
tive for  2 > 0 and  v > -1. 

Thus,  Equation  (2)  may  be  solved  by  transforming  it  into  a Riccati 
equation  and  then  into  a modified  Bessel  equation.  Once  boundary 
conditions  are  stipulated,  the  integration  constants  can  be  determined 
and  the  solution  for  u can  be  transformed  back  to  F.  This  is  not 
carried  out  for  the  general  case,  since  the  order  of  the  modified  Bessel 
functions  (v)  depends  on  the  value  of  C2.  I.e.,  it  is  preferable  to 

know  whether  v is  an  integer  or  non- integer  before  proceeding  with  the 
inverse  transformations. 

)lutions  of  this  type  are  applicable  to  equations  such  as  those 
steady  state  boundary  layer,  where  similarity  transforms  can  yield 
le  solutions. 


IV.  SOLUTION  FOR  F 

Utilizing  the  boundary  condition  that  when  the  radial  coordinate  is 
zero,  the  radial  velocity  is  zero  and  the  axial  velocity  represents  some 
mean  value  implies^ 


at  n = 0;  F = 0 and  f'  =0, 
or  from  Equation  (5),  C2  = 0. 

In  this  case  v = ±1,  Ij^(z)  and  I_i(z)  are  linearly  dependent  (in 
fact,  equal)  and  the  general  solution  (12)  becomes 

u = AIj^(z)  + BKj(z).  (13) 

Then , 

u{C)  = * 


^Handbook  of  Mathematiaal  Functions  With  Formulas,  Graphs ^ and  Mathe- 
matiaal  Tables^  Edited  by  M.  Abvamowitz  and  I.  Stecfun,  National  Bureau 
of  Standards^  Applied  Mathematias  Sevies,  Sf>j  June^  1964,  p.  374. 
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and  finally. 


F = 2 n 


results.  Note  that  the  numerator  and  denominator  of  the  expression  in 
brackets  may  be  divided  by  A or  B,  so  A and  B are  not  independent.  One 
constant  can  always  be  eliminated.  This  results  from  the  use  of  the 
Riccati  transformation  used  to  obtain  Equation  (10).  This  transforma- 
tion converted  a first  order  Equation,  (9),  into  a second  order  Equation, 
(10),  thereby  introducing  an  additional  constant.  The  two  constants  of 
integration  associated  with  Equation  (2)  are  Cj  and  A or  B.  The 
constant  Cj  appears  in  Equation  (2)  and  C2  has  already  been  evaluated. 

Finally,  the  expression  for  F may  be  written  as 


From  the  expression  for  F,  the  radial  and  axial  velocities  and  their 
derivatives  may  be  obtained,  in  a manner  analogous  to  Schlichting ' s 
treatment.  For  this  general  formulation,  boundary  conditions  must  be 
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specified  to  determine  Cj  and  A or  B. 
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